Emergent Collective Modes and Kinks in Electronic Dispersions 
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Recently, it was shown that strongly correlated metallic fermionic systems [Nature Phys. 3, 
168 (2007)] generically display kinks in the dispersion of single fermions without the coupling to 
collective modes. Here we provide compelling evidence that the physical origin of these kinks are 
emerging internal collective modes of the fermionic systems. In the Hubbard model under study 
these modes are identified to be spin fluctuations which are the precursors of the spin excitations 
in the insulating phase. In spite of their damping the emergent modes give rise to signatures very 
similar to features of models including coupling to external modes. 
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The description of nascent collective modes which 
emerge from elementary excitations on varying a control 
parameter g is an intensely studied field of research. The 
difficulty relies in the fact that in one limit of g the ele- 
mentary excitations dominate while in the other limit the 
collective modes dominate. In the vicinity of the transi- 
tion or around the crossover necessarily both degrees of 
freedom need to be taken into account so that the inter- 
play of both kinds of excitations is crucial. No simple 
theory assesses this interplay. 

Here we will focus on strongly correlated electronic sys- 
tems and especially on the metal-insulator transition in- 
duced by a repulsive interaction [/ on a lattice with a 
commensurate number of electrons per site. The sim- 
plest case is a local interaction with one electron per site 
on average [1]. For low values of U the electrons move 
through the lattice so that the system is metallic. For 
large values of U the hopping is blocked and the system 
is insulating with frozen charge degree of freedom. But 
the spin dynamics is still active. In leading order in t/U 
(t the hopping matrix element) this dynamics is captured 
by a Heisenberg model [2] . The collective modes are the 
spin excitations built from bound electron- hole pairs. 

When the system is still metallic, but close to its in- 
sulating regime, we intend to understand how the emer- 
gent spin modes influence the electronic quasiparticles. 
This issue is important to many strongly correlated sys- 
tems. One prominent example is high temperature su- 
perconductiviy where a large number of theories explains 
the attractive interaction between charge carriers by the 
interplay with spin fluctuations. One line of argument 
links the kinks that are observed in the dispersion of the 
fermionic holes, see for instance [3, 4, 5, 6, 7], to the in- 
teraction with bosonic modes. This is the usual reason- 
ing for phonons coupled to electrons [8]. Other bosonic 
modes, however, will engender the same sort of kinks, 
for instance plasmons [9]. In the high-Tc materials, spin 



fluctuations have an important influence on the quasipar- 
ticles, see e.g. Ref. [10]. They are likely candidates for 
the bosonic modes, see e.g. Ref. [11] where this is worked 
out in the fluctuation-exchange approximation. 

Byczuk et al. [12] recently showed by a sophisticated 
analysis of the equations of dynamic mean-fleld theory 
(DMFT) [1] that kinks in the electronic dispersion are a 
generic feature of strongly correlated electronic systems 
where the repulsive interaction is of similar strength as 
the kinetic energy. They stress that no coupling to a 
bosonic mode is needed. Indeed the model they study 
does not comprise any explicit bosonic mode; it is a 
fermionic Hubbard model. For particle- hole symmetric 
models dominated by the local self-energy the position 
of the kink was related by Byczuk et al. [12] to the quasi- 
particle weight 
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In the present work it is our aim to elucidate the phys- 
ical origin of the kinks. We provide evidence that the 
kinks result from the coupling to the bosonic resonance 
which is the precursor of the spin modes in the insulator. 
Thus we conjecture that the kinks in strongly correlated 
fermionic systems are induced by coupling to internal 
bosonic modes. The signature is very similar to the cou- 
pling to external bosons such as phonons [8] . Our finding 
also sheds light on signatures of spin modes in the elec- 
tronic dispersions of high-temperature superconductors. 

Our computation is also based on DMFT. This ap- 
proach reduces the extensive lattice problem to a self- 
consistency problem involving a single-impurity Ander- 
son model (SIAM) [1]. The latter can be viewed as 
an interacting site coupled to a semi-infinite chain of 
non-interacting fermions [13, 14] which is solved by dy- 
namic density-matrix renormalization (D-DMRG) [15, 
16]. This combination of D-DMRG and DMFT repre- 
sents a powerful tool for investigating the T = one- 



FIG. 1: (Color online) Panel (a) illustrates the relation be- 
tween real and imaginary part for an analytic ansatz of the 
real part with a kink: ReS = ^ (2^+l^-<»H.^+^<.|)(^"-i>") 
with A ^ 0.62177D, a = 0.15D, b = 0.7D; the imaginary part 
is computed by the Kramers-Kronig relation. Panel (b) shows 
the real and imaginary part of T,{uj) in DMFT at U = 2.0D. 

particle propagators of interacting lattice models [17-20] 
Its particular merit is to have a well-controlled energy 
resolution over the whole energy range [15, 21]. 

The model under study is the simplest displaying an 
interaction driven metal-insulator transition, namely the 
half-filled Hubbard model 

^ = E 4.S;. + ^E("'';T-l/2)(n,:;i-l/2). (2) 

At low values of U the ground state is metallic; above 
Uc2 ~ 3Z) the insulating phase becomes the ground state 
[1, 19, 22, 23]. 

Our analysis is facilitated by the direct numerical cal- 
culation of the local proper self-energy 2(0-"). This is done 
with the help of the improper self-energy 

Q{cj) ((d,(n_<, - l/2)|(n_, - 1/2)4)) (3) 

where we use the notation {{A\B)) for the Fourier trans- 
form of the time-dependent fermionic Green function 
—i{{A{t),B{0)}). If one considers doping the term 1/2 
in (3) is to be replaced by the average filling per site. 

Starting from the result S(w) = UF{lo)/G{uj) by 
Bulla et al. [24] we apply the Liouville operator in the 
equations of motion once more [25] yielding F{uj) = 
UQiLu)Go{uj) wherein F{lo) {{d„in_^ - 1/2)|4)). 
Substituting F{lj) by UQ{Lj)Gn{uj) and expressing G{lj) 
by Dyson's equation G~^{u!) = G~^{u!) — S(ijj) yields 

SH = c/2gH/(i + [/2Q(cj)GoM). (4) 

This expression is advantageous to use for small to mod- 
erate values of [/ ~ 2Z) where the computation of S(ijj) 
from the difference between the inverse bare and full 
propagators is numerically not reliable [20] . 

Figs. 1 and 2 display the generic behavior found for 
U not too far from the metal-insulator transition. Fig. 
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FIG. 2: (Color online) — ImE(cj) on a larger scale with two 
Fermi-liquid fits. The shaded region illustrates additional de- 
cay indicated by the arrows. 

1 (a) shows that a kink in the real part of a self-energy is 
linked to a trough-like feature in the imaginary part. This 
is a purely mathematical fact stemming from Kramers- 
Kronig relation. Panel (b) depicts real and imaginary 
part for a realistic self-energy as it results from the 
DMFT calculation. The kinks in the real part and 
the trough in the imaginary part are clearly discernible 
though not as neatly as in the analytic function of panel 
(a). This comes from small spurious wiggles in — ImE re- 
sulting inevitably from the dcconvolution of the DMRG 
raw data [26]. 

In Fig. 2 we address the physical meaning of the 
trough-like feature. There are two ways to understand 
it based on Fermi liquid theory. 

(i) The trough itself, ranging approximately from 
— O.ID to OAD, is fitted by a narrow curve (dashed- 
dotted fine). Outside the trough — ImE(a;) is then much 
lower than the extrapolated fit. Since — ImE(w) is the 
decay rate for the quasiparticlcs, this would be much 
lower relative to its extrapolated value. We do not see 
a good reason for such a decrease of the decay because 
the decrease due to reduced phase space for three quasi- 
particles should occur beyond about three times ZD/2 
which is sa OAD [19, 20, 22], i.e., significantly larger than 
the extension of the trough. 

(ii) The Fermi liquid theory extends to higher values, 
for instance « 0.5D as is still consistent with the above 
crude estimate, so that fits such as the dashed one in Fig. 

2 are justified. Indeed, the fit works very nicely with a 
moderate coefficient for the quartic term. This view im- 
plies that around 0.15D additional decay becomes pos- 
sible which extends up to OAD. An additional decay 
channel is well possible. It sets in only above a certain 
energy because excitations of a certain minimum energy 
are involved. 

So, among the two hypotheses we favor the second one. 
It explains also rather naturally why the quadratic coef- 
ficient is so low in spite of the very narrow trough. 
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FIG. 3: (Color online) Deconvolved imaginary part of the 
local spin susceptibility at positive frequencies for various in- 
teractions U in the metallic phase. 



We are aware that from a puristic point of view on 
Fermi liquid theory its apphcabihty ends at the borders 
of the narrow trough as described above in (i). We do 
not claim that this view is invalid. But we advocate the 
alternative view (ii) because it provides an intuitive way 
to understand the self-energy behavior at low energies 
in terms of quasi-particles coupled to emergent collective 
modes. This coupling is the origin for the deviations from 
the dashed curve in Fig. 2. The latter is regarded as an 
effect on top of the underlying Fermi liquid description 
resulting from the additional decay channel. This is in- 
duced by scattering from an emergent collective mode 
which has to be identified. 

Given the fact that it becomes important only for fi- 
nite, though small energies we aim for a mode which is 
dominated by such a finite, though small energy. Fur- 
thermore, it exists only close to the metal-insulator tran- 
sition. We shall see that its energy decreases towards the 
transition U —f Uc2- Since the insulator is a paramag- 
net with disordered local spin moments [1, 27] a natural 
candidate are the emergent spin fiuctuations. 

In the framework of the limit of infinite dimension 
d ^ oo the propagation of a collective mode from site 
i to site j scales like d^'*^^' where | • | stands for the 
taxi cab metric. Hence the collective modes are almost 
dispersionless and thus local. Only for particular wave 
vectors which are of measure zero a non-local propaga- 
tion makes itself felt, see, e.g., Ref. 28. In the complex 
diagrams describing the single particle motion the prop- 
agation of collective modes (particle-hole pairs) occurs in 
such a way that it is summed over. No particular mo- 
menta of measure zero matter. Thus it is fully sufficient 
to investigate the local response. 

In Fig. 3 the local spin susceptibility Xspinl'^) is shown 
which we have computed for positive frequencies denot- 
ing it by x^in('^)- It is obtained from the effective 
SIAM occurring in the self-consistency loop of the DMFT 
[1, 20]. In the SIAM it is determined as the local sus- 
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FIG. 4: (Color online) Kink positions o^kink as derived from 
the quasi-particle weight Z via Eq. (1); Z is found either 
from the propagator G [19] or from the self-energy E(aj) in 
Z = (l-ac^E(0))~\ Most directly, a fit Aoj + B{\u} -LO^nkl - 
|cli + ti^kinkl) to ReE(Lj) is used for tJkink; 'i'kink IS Compared to 
the energies where x^pin shows a peak at low \ll>\. 



ceptibility at the head of the chain. Numerically we 
employ again D-DMRG for some broadening j] which is 
then eliminated by deconvolution [26]. This deconvolu- 
tion gives rise to some uncertainty in the shape of the 
frequency dependence of the susceptibility. 

The excitation operator is 25"^ ^ n^- — at the chain 
head. A strongly pronounced peak catches the eye. Its 
peak energy moves towards to — for U —s- Uc2- In paral- 
lel, its height increases such that its total weight tends to 
a finite value [27] . This peak is the precursor of a (5 peak 
at zero energy in the paramagnetic insulator. There it 
reflects the fact that a spin can be rotated without any 
cost of energy. Still in the metallic phase, the peak is 
a resonance made from an almost bound quasiparticle 
and a hole. It has some width because it may decay into 
scattering states of its constituents. 

From the data for Imx^jn we deduce the peak position 
by fits assuming two lorentzians to account for the asym- 
metry of the peak shape. The lorentzians are multiplied 
by factors tanh(a;/ti>o) to account for the linear vanishing 
of Imx^in(ti') for uj ^ 0. The relevant peak position is 
the one of the lorentzian with more weight. The error 
bars account for the uncertainties related to the details 
of the fit procedure, e.g., for U « 1.9D where the weight 
appears to be distributed equally over both lorentzians. 

The results are compared in Fig. 4 with the kink posi- 
tions which were determined in several ways. We use 
Eq. (1) to deduce the kink position from Z which in 
turn is determined either by Z = (1 — 9„S(0))~^ or by 
Z-^ = D'^d^G{0)/2 [19]. Or the kink position is deter- 
mined directly by a fit to Rel](0) (see caption of Fig. 4). 
The three ways to find the kink positions agree very well 
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providing consistent data. 

The peak positions agree remarkably well with the kink 
positions. In particular for large values of U the agree- 
ment is striking. It is for these larger values U ^ 2D 
that both the kink and the peak in the susceptibility are 
clearly discernible and well-defined. So we deduce that 
the additional decay channel seen in Fig. 2 results from 
the excitation of the spin resonance by the propagating 
single fermionic quasiparticle. This finding strongly sup- 
ports our claim that the kink is in fact due to emergent 
internal modes. Here these modes are the spin fluctua- 
tions which develop already in the metallic phase. 

Thereby, an intuitive physical picture of the origin of 
the kinks is found. One major advantage of this picture 
is that one can transfer it to finite dimensions where the 
collective modes are dispersive so that the momentum 
dependence matters. Kinks are to be expected where 
momentum and energy conservation for the scattering of 
a quasiparticle from a collective mode is fulfilled. 

In conclusion, we have provided compelling evidence 
for a link between the kinks in fermionic dispersions in 
strongly correlated systems and emergent internal col- 
lective modes, spin fluctuations in particular. We agree 
completely with the phenomenon established by Byczuk 
et al. [12]. But our physical picture of the phenomenon is 
different because we view the kinks as the consequence of 
inherent bosonic modes. An important concomitant as- 
pect is that the Fermi liquid theory does not break down 
already at the scale of Wkink- It extends to about 2ZD 
where Z is the quasiparticle weight. 

Our finding provides important information on the 
possible interpretation of kinks in electronic dispersions 
in many strongly correlated systems and in cuprate sys- 
tems as they occur in high-temperature superconductors 
in particular. Such kinks can be the consequence of 
emerging bosonic modes, i.e., resonances even if these are 
still strongly damped. For instance, qualitative support 
is provided to results based on the fluctuation-exchange 
approximation for cupratcs [11] Moreover, the coupling 
between the single particles and the collective modes is 
generically substantial. Certainly, further investigations, 
for instance away from half-filling, are called for. 
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